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in ■ Abstract. As so far, the redshift of Gamma-ray bursts (GRBs) can extend to z ~ 8 which 

^^ I makes it as a complementary probe of dark energy to supernova la (SN la). However, the 

calibration of GRBs is still a big challenge when they are used to constrain cosmological 
models. Though, the absolute magnitude of GRBs is still unknown, the slopes of GRBs 

O ■ correlations can be used as a useful constraint to dark energy in a completely cosmologi- 

cal model independent way. In this paper, we follow Wang's model-independent distance 
measurement method and calculate their values by using 109 GRBs events via the so-called 

k> . Amati relation. Then, we use the obtained model-independent distances to constrain ACDM 

$H I model as an example. 

^ ! 
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1 Introduction 

Since the discovery of an accelerated expansion of our Universe through the observations of 
supernova la (SN la) [1, 2], many cosmic observations have been used to explore the nature of 
dark energy which has negative pressure and pushes the Universe into an accelerated phase. 
Particularly, the redshift of Gamma-ray burst (GRBs) can extend to higher redshift z ~ 8. 
This merit makes it as a complementary cosmic probe to SN la. However, a big challenge, 
the so-called circular problem, has to be overcome when one calibrates the GRBs to obtain 
the distances at different redshifts. In [3], Schaefer et. al. have derived the distance modulus 
in ACDM model by using five GRBs correlations. The so-called circular problem would 
be committed when one uses the resulted distance redshift relation to constrain any other 
cosmological models beyond ACDM. To overcome this circular problem, Li, et. al [4] put the 
GRBs correlation and cosmological model constraint together. Via the Markov Chain Monte 
Carlo (MCMC) method, they fixed the calibration parameters and constrained the model 
parameter space simultaneously. However, it looks like using a 'ruler' having no marks to 
measure the length of an object. And after the measurement, the length and scale of ruler 
are given together. So the lack of calibration makes the GRBs lose the power to constrain 
cosmological models. Cosmography method was considered in [5] by parameterizing the 
luminosity distance d^ in terms of deceleration qq, jerk jq and snap sq parameters. In this 
way, the cosmological model dependent problem was removed. Liang et. al. [6] calibrated 
GRBs by using low redshift SN la and obtained a tight constrain to the cosmological model 
parameter space. This method was reconsidered by Wei [7, 8]. By analyzing the calibration 
process carefully, one can find the potential drawback. The calibrated GRBs correlation 
makes GRBs have the same luminosity distance-redshifts relation as SN la at lower redshifts. 
In this way, the relation is extended to higher redshifts. Equivalently, one just extends the 
luminosity distance-redshifts relation of SN la to higher redshifts. So, it makes the obtained 
luminosity distance-redshifts relation strongly depends on that of SN la. The worst thing 
is that it makes the data points of GRBs useless, because one has known the luminosity 
distance-redshifts relation at high redshifts from SN la, though no SN la is found at the 
high redshift regions of GRBs. Alternatively, Wang presented a model-independent distance 
measurement from GRBs calibrated internally [18]. The main point of Wang's method is 



that the statistical errors of correlation parameters aa, o"fe and systematic error agys obtained 
in ACDM models are used, but are not the correlation parameters a and b themselves, for 
the definitions of a, h please see Eq. (2.2). The viability of this implement comes from the 
observations that the errors of correlation parameters are almost the same for different values 
of VtmQ for ACDM model, though the values of a and h are really different. Then, in terms of 
a set of model-independent distance measurements, the cosmic constraint from GRBs is set 
up via cubic spline interpolation from cosmological model independent distance ratio fp{zi). 
The merits of this method are follows: (i) the constraint from GRBs is in a cosmological 
model independent way. It alleviates the circular problem. So, it can be used to constrain 
any other cosmological models, (ii) It is not calibrated by any other external data sets. It 
does not suffer any consistent problem when it is combined with other data sets as cosmic 
constraints, (iii) The cosmological model independent calibration is done firstly. It means 
that the 'ruler' has been marked, (iv) Though the absolute magnitude of GRBs is unkown, 
the slopes of GRBs correlations can be used as cosmological constraints. 

Recently, Wei [8] used 109 GRBs data points via Amati relation [11-13] calibrated by 
SN la to constrain cosmological model. Based on the points mentioned above, to alleviate the 
data sets dependence and circular problem, we shall present our calculation results based on 
Wang's method via Amati correlation. In fact, we find a new Xgrb which only depends on 
the slopes of GRBs correlation and make the absolute magnitude of GRBs irrelevant. For the 
details, please see section 2. As results, a set of model-independent distance measurements 
are obtained on the basis of Amati relation that can be used to constrain cosmological models. 

However we have to stress that the Amati relation has been criticized for many reasons 
in the literatures. Li demonstrated an ambiguity in determining the redshifts of GRBs [14]. 
Of course, this ambiguity can be overcome when the redshift is well determined. The major 
criticism came from Nakar and Prian [15] who developed a test for the Amati relation even in 
the case where the redshifts of GRBs were unknown. The test was also generalized by Band 
and Preece [16]. They concluded that the Amati relation suffered the problem of selection 
bias. Recently, the authors of [17] have also concluded that the Amati relation is an artifact 
of selection effects with the burst population and the detector. They also point out that 
the Amati relation is failed whether or not the bursts have measured spectroscopic redshifts. 
If this is true for Amati relation, the results obtained based on Amati relation would be 
unreliable. But it is still in debating. So, before the dust settles down, we still assume that 
Amati relation is reliable in this paper. 

This paper is structured as follows. In section 2, the values of cosmological parameter 
dependance in ACDM model are presented as shown in [18]. The errors of correlation pa- 
rameters (Ta, (Tfo and systematic error a^yg are also calculated via Amati relation with 109 
GRBs data points. The 5 bins model-independent distance measurements will also be found 
in this section. In section III, we use the resulted data points to constrain ACDM model as 
an example. A summary and discussion are put in section IV. 

2 Calibration of GRBs and Model-independent Distance Measurement 

Following the work of [3], we consider the well-known Amati Ep^i — Eiso correlation [9, 11-13] 
in GRBs, where Ep^i = Ep^obsO- + ^) i^ the cosmological rest-frame spectral peak energy, and 
Eiso is the isotropic energy 

Eiso = ^TidlSi,oio/{l + z) (2.1) 



in which di and Sboio are the himinosity distance and the bolometric fluence of the GRBs 
respectively. Following [3], we rewrite the Amati relation in the form of 

log^ = a + 61og-|^. (2.2) 

erg SOOkeV 

One fitts the Amati relation through the minimization of x^ which is given by [3] 
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x^ = log _3b^ (2.4) 
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where di is defined as [18] 

dL = Ho{l + z)r{z)/c, (2.6) 

and the errors are calculated by using the error propagation law [10]: 
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a^i= , Z (2.7) 

lnlO£;p,i ^ ' 



'^Sl^olo.; 



In li)Shoio-i 



{21 



By calculating the value of y^, we find it is large and dominated by the systematic 
errors, and on the contrast the statistical errors on a and h are small. Following [3], the 
systematic error Ugys can be derived by required \^ = v (the degrees of freedom). At last, 
the total error a■'^^^ = (7%^^ + alyg is obtained. It would be noticed that in our case, the 
best fit value of a will be less 21og(c/i?o) than that in the definition of luminosity distance 
di = 0- + z)r{z) [18]. For this definition, the value of Hq is absorbed into the calibration of 
GRBs because of the lack of fixing the absolute magnitude of GRBs. Then, in our treatment 
the results will be Hq free. 

As shown in [18], the calibration of GRBs is cosmological model dependent, because the 
values of a and b are obtained on the basis of ACDM model. So the calibrated result can not 
be used to constrain any other cosmological models. This the so-called well-known model 
dependent problem or circular problem. In Tab. 1, the cosmological parameter dependence 
are shown. With a careful observation, one can find that the Icr errors bars of o, b and system 
are almost the same which do not depend on the cosmological parameters. In this paper as 
suggested by Wang [18], we shall only use the systematic error and the errors of a and b in 
the case of OmO = 0.27 for ACDM model as the standard values. That is to say the values of 
(Ta, (Jb and Gsys will be used in the following sections not the values of a and b derived from 
ACDM model. 

Following the work of Wang [18], the Xgrb '^^ ^ cosmological model is given by 

Xgrb - /^ rjTZZlTTd;!^^. ' ^'^■^1 
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0.27 



-3.392 ± 0.0368 

1.583 ± 0.0729 

0.324 



a, 
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0.2 



-3.343 ± 0.0358 

1.600 ±0.0744 

0.328 



a 



mO 



0.4 



-3.467 ±0.0364 

1.554 ±0.0725 

0.321 



Table 1. Systematic error and values of a and b for GRBs Amati relation in the cases of ilmo 
0.27, 0.2, 0.4 in ACDM model. 
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To constrain a cosmological model, one uses a set of model-independent distance mea- 
surements {fp(zi)}: 






(2.13) 



where r(z) = dL{z)/{l + z) is the comoving distance at redshift z, zq is the lowest GRBs 
redshift^. Here, the definition of r^ is different from Wang's definition [18] rp{z) = (1 ± 
z)^'"^ HQr{z) / [zch) where h = i:fo/(100kms~ Mpc~ ). In our definition, the distance mea- 
surement fp{zi) and rp{z) is completely Hq free. It can be seen from the facts that the 
definition of r{z) is r{z) = c/Hq J^ dz' /E{z') where H'^{z) = HqE'^{z) is the Hubble param- 
eter. In terms of our definition, d^ can be rewritten as 



di = z{l + z)^/^rp(zo)rp(z). 



(2.14) 



We divide the redshifts of GRBs into N bins, i.e. {zi},i = 1, ..., A^, and assume the corre- 
sponding values of {rp{zi)},i = 1,...,N which do not depend on any cosmological models. 
Then the values of rp(z) at arbitrary redshift z can be obtained by cubic spline interpola- 
tion from {fp{zi)}. So, the values of ^^(z) and diiz), etc at redshift z can be found easily. 
Given each set of {rp(zj)},z = 1, ■■■,N, we calibrate the GRBs and calculate the likelihood 
simultaneously via Eq. (2.9) by using Markov Chain Monte Carlo (MCMC) method [19]. 
The MCMC is a global fitting method which is used to determine the cosmological param- 
eters. In adopting the MCMC approach, we generate using Monte Carlo method a chain 
of sample points distributed in the parameter space according to the posterior probability, 
using the Metropolis-Hastings algorithm with uniform prior probability distribution. In the 
parameter space formed by the constraint cosmological parameters, a random set of initial 
values of the model parameters is chosen to calculate the x^ or the likelihood. Whether 
the set of parameters can be accepted as an effective Markov chain or not is determined by 
the Metropolis-Hastings algorithm. The accepted set not only forms a Markov chain, but 
also provides a starting point for the next process. We then repeat this process until the 
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0.17 was used in [18]. In this work, the lowest redshift of GRBs is zq = 0.0331 



established convergence accuracy can be satisfied. The convergence is tested by checking the 
so-called worst e- values [the variance(mean)/niean(variance) of 1/2 chains] R—1 < 0.005 [19]. 
As results, we obtain a set of distances {fp{zi)},i = 1, ...,N which are independent on any 
assumption of cosmological parameters. It comes from the observations that, in the process 
of calibration, the statistical and systematic errors of cja, a^ and agys are only used. We do 
not use any values of a and b calibrated in ACDM model. The important thing is that the 
statistical and systematic errors aa, Cf, and dsys are almost model parameter independent. 
Thanks to this feature, this method is model-independent. Here, in the MCMC analysis, 
we take a, b and A^ {rp{zi)},i = 1, ..., A as free parameters. So, the degree of freedoms is 
v = 109 — 2 — N. Once these values of {fp{zi)},i = 1, ..., A are obtained, a cosmological 
model can be constrained by GRBs via the x^ 

Xgrb = [^rp{zi)] ■ {Cov^]^^)ij ■ [Afp{z,)], (2.15) 

Arp{z,) = f^p''"'izi)-rpiz,), (2.16) 

where fp{zi) is defined by Eq. (2.13) and {CovQj^^)ij,i,j = 1...A is the covariance matrix. 
In this way, the constraints from larger observational GRBs data points are projected into 
relative smaller number of points. Of course, this method can be generalized to discuss other 
problems. 

Now, we present some discussion about the treatment of rp{zo)- Here, we do not calcu- 
late the value of rp at the redshift zq via its definition (2.13) according to any cosmological 
model. Because, if we calculate the value for any cosmological model, a model or cosmolog- 
ical parameter dependence will be introduced again. The values of rp{z()) are fixed by the 
calibration relation (2.2), i.e. via the relation 

rAzo) = W'^/'i^^f' ,, ^\ ,,,,, . (2.17) 
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In fact, once this relation is used, one can recast Eq. (2.9) into the following form 

NgRB /y _ Yf 

Xgrb = 2_/ ('data\2 ' (^-^^^ 

where Yi is defined as 

Yi = blogEp^i{zi) - log Sboioizi) - 21ogZi - 2log fp{zi) (2.19) 

and Yq = Yi{z(j), the Cy"/" is the total la errors of data sets 



One can see that the parameter a is removed from this new form of Xgrb- ^^ other words, the 
information about the slope b of GRBs correlation is used alone. So the absolute magnitude 
a of GRBs is irrelevant in this method. Of cause, the rp{zQ) can be fixed by consulting a 
special cosmological model or other data sets. But in that way, the circular problem and 
data sets relevance problem will come back. Equivalently, the absolute magnitude of GRBs 
was fixed in a special cosmological model. So, for every possible values of b in every running 
of MCMC, the Xgrb ^^ calculated. Here, we must keep in mind that the values of parameter 



b are not taken from Tab. 1. It is treated as a free parameter. So, this method is completely 
cosmological parameter in-dependent and self consistent. 

We divide the redshifts into N = 5 bins and run the MCMC codes which is based on the 
publicly available CosmoMC package [19]]. The chains have worst e- value R — 1 = 0.0017 
which is much smaller than 0.005. The resulted model-independent distances and covariance 
matrix from 109 GRBs are shown in Tab. 2 and Eq. (2.22). As already mentioned above, 
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Table 2. Distances measured form 109 GRBs via Amati relation with Ict upper and lower uncertain- 
ties. 

zq = 0.0331 is adopted in this work. The {fp{zi)},i = 1, ...,5 correlation matrix is given by 

/ 1.0000 0.7780 0.8095 0.6777 0.4661 \ 

0.7780 1.0000 0.7260 0.6712 0.3880 

(Covgrb) = 0.8095 0.7260 1.0000 0.6046 0.5032 , (2.21) 

0.6777 0.6712 0.6046 1.0000 0.1557 

V 0.4661 0.3880 0.5032 0.1557 1.0000/ 



and the corresponding covariance matrix is given by 



where 



{CovGRB)ij = (y{rp{zi))a{rp{zj)){CovGRB)ij, 



a{fp{z.,))=a{fp{z,))+, if fp{z)>f.,{zf-'-- 
a{rp{zi)) = a{fj,{z{))- , if f.p{z) < fp{zf^'\ 



(2.22) 
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here the a{fp{zi))'^ and (j{fp{zi))~ are the la errors listed in Tab. 2. The marginalized 
values and the corresponding upper and lower bounds are used when GRBs is used as a 
cosmological constraint. The ID distributions of model parameters are shown in Fig.l. 
The distance measurements from 109 GRBs via Amati correlation with \a error bars are 
shown in Fig. 2, where the solid lines correspond to ACDM model with different values of 
^mo = 0.2, 0.27, 0.4 from up to bottom respectively. 

The data points shown in Tab. 2 and Fig. 2 are totally cosmological model independent, 
so one can use them to constrain any other cosmological models. In the next section, as an 
example, these obtained data points will be used to constrain ACDM model. From Fig. 2, 
one can find that GRBs favors large values of VtmO for ACDM model. This point can be 
confirmed in the next section. And the concordance model {Q-mO = 0.27) is almost at the 
boundary of la regions of data points. So the null hypothesis that Amati relation is based 
on cosmology can be rejected slightly greater than 4.5it as shown in Fig. 2. 
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Figure 1. The ID marginalized probabilities of model parameters. 




Figure 2. The mean values of distance measurement from 109 GRBs via Amati relation with la 
error bars. The lines from up to bottom correspond to ilmo = 0.2,0.27,0.4 respectively. 



3 Cosmological Constraint to Dark Energy 

In our calculations, we have taken the total likelihood function L oc e~^ ' ^ to be the products 
of the separate hkelihoods of SN, BAO, CMB and GRBs. Then we get x^ 



X — XsN + Xbao~^ Xcmb + Xgbbi 



(3.1) 



where the separate likelihoods of SN, BAO, CMB are shown in the Appendix A. The Xgrb 
is the form of Eq. (2.15). In this work, we only consider the ACDM model as a simple 
example. Its generalization to constrain other cosmological models is straight forward. In 
the case, when SN la and GRBs are combined as cosmic constraint, we have Xmin ~ 547.727 
and the best fit values of model parameter $7^0 = 0-274tooi94- If the SN la is used alone, 

o _l_0 OOOO 

one has Xmin ~ 542.680 and ilmo = 0-270_q'q2]^3. The corresponding contour plot is shown 
in Fig. 3. When combing the cosmic observations considered in this work together, we 
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Figure 3. Color Online: 2D contours plots in the 51a ~ ^m plane. The solid blue lines denote the 
1(7 and 2(7 regions from SN Union 2 alone. The red dashed lines denote the Itr and 2a regions from 
the combination of SN Union 2 and 109 GRBs. 



have the resulted x. 
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549.383 and n 



mO 



0.279l[]:|5}^^. Also, the best fit values of 



Qrno = 0.277J:o:o}28 ^nd Xmin = 544.451 when GRBs is not used for comparison. The 2D 
contour plots are shown in Fig. 4. We can find that when GRBs is used, the relative 
errors of model parameters is shrunken. Here, we can compare our result with Wei's one 
where he used 50 low redshift GRBs to calibrate the Amati relation. So, we add 59 high 
redshift data points of GRBs to SN Union 2 data sets. After constraint via MCMC, we have 
Xmin — 565.919 and QmO = 0-271^ooi92- "^^^ corresponding 2D contour plot is shown in Fig. 
5. The constrained results show that relative large values of O^o are favored when GRBs 
data points are employed. This is consistent with the clues shown in Fig. 2. 
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Figure 4. Color Online: 2D contours plots in the JIa — ^m plane. The solid blue lines denote the la 
and 2(T regions from SN+BAO+CMB. The red dashed lines denote the la and 2a regions from the 
combination of SN+BAO+CMB+GRBs. 



4 Summary and Discussion 

In this paper, by using 109 GRBs data points via Amati relation, we have derived five data 
points of distance measurements which do not depend on any cosmological models, i.e. in a 
model-independent way, based on the method firstly advocated by Wang [18]. Then it can be 
used to constrain any other cosmological models without the so-called circular and data sets 
consistence problem. We also find out that the concordance model (0^0 = 0.27) is almost at 
the boundary of la regions of data points. So, one can reject the null hypothesis that Amati 
relation is based on cosmology slightly greater than A. 5a as shown in Fig. 2. 

When GRBs (the five data points of distance measurements) is used as a complemen- 
tary cosmic constraint to ACDM model as an example, the errors of model parameters are 
shrunken. Though the constraint is not much tighter than that calibrated via low redshift 
SN, the difference is very small. And, the results can be used to constrain other cosmological 
model beyond ACDM. The important thing is that it does not depend on SN la data set. 
The merits of this method has been mentioned in the introduction. Here, we list possible 
potential drawbacks of this completely model in-dependent method as follows: (i). The con- 
straint is not much tighter than that obtained from calibration by using SN la data points. 
The looser may come from the lack of enough GRBs data points or the number of bins is 
not larger, (ii). How to design a cosmological model indicator, here it is model independent 
distance measurement fp{z), which is sensitive to distinguish the cosmological models may 
be the main points of this method. In this way, the lack of having no enough GRBs data 
points can be avoided in some senses, (iii) Another negative effect comes from the choice of 
fiducial redshift zq. Some average values of rp{zo) would be better. 
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Figure 5. Color Online: ID contours plots in the JIa — ^m plane. The solid pink lines denote the 
\a and 2a regions from SN Union 2 + 59 GRBs from Wei's results. The red dashed lines denote the 
ItT and 2(7 regions from the combination of SN+GRBs in our case. 

At last we have to alert the readers that our analysis is based on the assumption that 
Amati relation is reliable. However, in the literatures, the Amati relation has been criticized 
for many reasons as mentioned in the introduction. If that is true, our analysis would be 
unreliable. 
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A Cosmic Observations: SN, BAO and CMB 

A.l Type la Supernovae constraints 

Recently, SCP (Supernova Cosmology Project) collaboration released their Union2 dataset 
which consists of 557 SN la [21]. The distance modulus /x(z) is defined as 



ii'thi.z) = 51ogio['^L(2;)] + fJ-O, 



(A.l) 



where dL{z) is the Hubble-free luminosity distance H^dL^z) / c = HQdA{z){l + z)"^ /c, with Hq 
the Hubble constant, defined through the re-normalized quantity h as Hq = lOO/i km s~ Mpc~ , 



10 



and /io = 42.38 — Slog^g h. Additionally, the observed distance moduli fJ-obsi^i) of SN la at 
Zi is 

fJ'obsizi) = niobsizi) - M, (A.2) 

where M is their absolute magnitudes. 

For the SN la dataset, the best fit values of the parameters ps can be determined by a 
likelihood analysis, based on the calculation of 

2/ ^A _ y^ {fJ-obsjZi) - fJ.th{Ps,Zi)}'^ 
SN ^i 

_ ^^ {5logio[dLiPs,Zi)] -mobs{zi) + M'} 

SN « 

where M' = /xq + Af is a nuisance parameter which includes the absolute magnitude and the 
parameter h. The nuisance parameter M' can be marginalized over analytically [22] as 



X (Ps) = -21n / exp 
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Relation (A. 3) has a minimum at the nuisance parameter value M' = B /C, which contains 
information of the values of h and M. Therefore, one can extract the values of h and M 
provided the knowledge of one of them. Finally, note that the expression 

xIn{Ps,B/C) = A-{B^/C), 

which coincides to Eq. (A. 4) up to a constant, is often used in the likelihood analysis [22, 23], 
and thus in this case the results will not be affected by a flat M' distribution. 

A.2 Baryon Acoustic Oscillation constraints 

The Baryon Acoustic Oscillations are detected in the clustering of the combined the 2dF 
Galaxy Redshift Survey (2dFGRS) and Sloan Digital Sky Survey (SDSS) main galaxy sam- 
ples, and measure the distance-redshift relation at z = 0.2. Additionally, Baryon Acoustic 
Oscillations in the clustering of the SDSS luminous red galaxies measure the distance-redshift 
relation at z = 0.35. The observed scale of the BAO calculated from these samples, as well 
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as from the combined sample, are jointly analyzed using estimates of the correlated errors to 
constrain the form of the distance measurement Dv{z) [24-26] 



z 



dz' ,\"3 



where sinn(x) = sin(x),x,sinh(x) for il^ < 0, $1^ = and [7^ > respectively. The peak 
positions of the BAO depend on the ratio of Dv{z) to the sound horizon size at the drag 
epoch (where baryons were released from photons) Zd, which can be obtained by using a 
fitting formula [27]: 

with 

hi = 0.313(O„/i2)-0-4i9[i + 0.607(a;„/i2)0-674]^ (A.7) 

62 = 0.238(O„/i2)0-223. (A.8) 

In this paper, we use the data of rs{zd)/Dv{z) extracted from the Sloan Digitial Sky Survey 
(SDSS) and the Two Degree Field Galaxy Redshift Survey (2dFGRS) [28], which are listed 
in Table 3, where rs{z) is the comoving sound horizon size 

, , /■* Csdt r Csda p , c, 

rs{z) =c / = c / ^— = c / dz 



a Jq a?H J^ H{z) 

V(i+-) da 



V3J0 a'^H{a)y^lTWhIWhWJ 

where Cs is the sound speed of the photon— baryon fluid [29-31]: 



(A.9) 



9 4 Ph(z) 4 , ilj, , , . s 



and here 0^ = 2.469 x 10~5/i"2 for Tcmb = 2.75K. 



z 


r,{zd)/Dv{z) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



Table 3. The observational rs{zd)/Dv{z) data [25]. 
Using the data of BAO in Table 3 and the inverse covariance matrix V^^ in [25]: 

1 _/ 30124.1 -17226.9 \ , . 

~ V -17226.9 86976.6 / ' ^ ' ^ 

Thus, the Xbao(Ps) ^^ given as 

xlAoiPs)=^'V-'X, (A.12) 
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where X is a column vector formed from the values of theory minus the corresponding 
observational data, with 

and X* denotes its transpose. 

A. 3 Cosmic Microwave Background constraints 

The CMB shift parameter R is provided by [32] 



Riz,) = ^/^5smn[^^ /" 4^], (A.14) 



here, the redshift z^ (the decoupling epoch of photons) is obtained using the fitting function 

[33] 

z, = 1048 [1 + omuAiribh'^y^-^^^] [i + gi{n,nh^y^] , 

where the functions gi and g2 read 

51 = 0.0783(J^b/i2)-0-238 (1 + 39.5(J7b/i2)0-763)-i , 
c,2 = 0.560 (1 + 21.1(fife/i2)i-8i)~^ . 

In additional, the acoustic scale is related to the first distance ratio and is expressed as 





7 - year ML 


7 — year mean 


error, a 


Ia{z,) 


302.09 


302.69 


0.76 


R{z.) 


1.725 


1.726 


0.018 


z^ 


1091.3 


1091.36 


0.91 



Table 4. The observational Ia: R, z* data [34]. The ML values are used in this work as recommended. 

Using the data of Ia,R,z^: in [34], which are listed in Table 4, and their covariance 
matrix of [Ia{z*),R{z^),z^] referring to [34]: 

/ 2.305 29.698 -1.333 \ 
C"^ = 29.698 6825.270 -113.180 , (A.16) 

V -1.333 -113.180 3.414 / 

we can calculate the likelihood L as Xcmb ~ — 21nL: 

xImb = Ad^[C-\d„dj)][Ad,]\ (A.17) 

where Adi = di — df"'^"' is a row vector, and di = {lA,R,z,f). 
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